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ABSTRACT 

In this work we analyse the properties of cosmic voids in standard and coupled dark energy 
cosmologies. Using large numerical simulations, we investigate the effects produced by the 
dark energy coupling on three statistics: the filling factor, the size distribution and the stacked 
profiles of cosmic voids. We find that the bias of the tracers of the density field used to identify 
the voids strongly influences the properties of the void catalogues, and, consequently, the pos¬ 
sibility of using the identified voids as a probe to distinguish coupled dark energy models from 
the standard ACDM cosmology. In fact, on one hand coupled dark energy models are char¬ 
acterised by an excess of large voids in the cold dark matter distribution as compared to the 
reference standard cosmology, due to their higher normalisation of linear perturbations at low 
redshifts. Specifically, these models present an excess of large voids with i?eff > 20,15,12 
h~^ Mpc , at z = 0, 0.55,1, respectively. On the other hand, we do not find any significant 
difference in the properties of the voids detected in the distribution of collapsed dark matter 
halos. These results imply that the tracer bias has a significant impact on the possibility of 
using cosmic void catalogues to probe cosmology. 
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1 INTRODUCTION 

Despite the fact that the presently accepted standard cosmo¬ 
logical model, the so-called ACDM scenario, appears to be 
fully consistent with most of the available observations (see e.g. 
Planck Collaboration et al. 2015a), it still presents some open is¬ 
sues in the detailed description of the distribution of matter at 
small scales. One of such properties that still appears problem¬ 
atic is the observed abundance of dwarf galaxies in the underdense 
regions of the Universe, which is found to be significantly lower 
than what predicted by large N-body simulations carried out within 
the ACDM cosmology. This problem, that was pointed out for the 
first time by Peebles (2001), goes under the name of the void phe¬ 
nomenon, and it has been discussed by several authors over the past 
years (see e.g. Tinker & Conroy 2009; Sutter et al. 2015a). 

Besides the poor theoretical understanding of a cosmo¬ 
logical constant as source of the observed accelerated expan¬ 
sion of the Universe (Weinberg 1989), the void phenomenon 
is therefore one of the few observational tensions that moti¬ 
vate the investigation of alternative cosmological scenarios, to¬ 
gether with the so-called cusp-core problem (de Blok 2010), the 
satellite problem (Bullock 2010), the too big to fail problem 
(Boylan-Kolchin, Bullock & Kaplinghat 2011), and the recently 


detected tension between the CMB- and cluster-based estimations 
of (78, the r.m.s. of the mass density field within a sphere of radius 
8 h~^ Mpc (Planck Collaboration et al. 2015b). 

A relevant class of alternative cosmological models that has 
been widely investigated in recent years is given by the so-called 
coupled dark energy scenario (cDE hereafter, see e.g. Wetterich 
1995; Amendola 2000, 2004; Farrar & Peebles 2004; Baldi 201 lb). 
In these models a dynamical scalar field sourcing the accelerated 
cosmic expansion (see e.g. Wetterich 1988; Ratra & Peebles 1988) 
is coupled to cold dark matter (CDM) particles resulting in a di¬ 
rect exchange of energy-momentum between these two cosmic 
components. Such interaction gives rise to a new fifth force act¬ 
ing on CDM particles, possibly capable to make the voids emptier 
(Nusser, Gubser & Peebles 2005). Other possible ways to address 
the void phenomenon have been proposed, such as, for example, 
a modification of gravity at very large scales (Li & Zhao 2009; 
Clampitt, Cai &. Li 2013; Spolyar, Sahlen & Silk 2013). 

The main effects of cDE models on the large-scale matter dis¬ 
tribution in the Universe, as well as on the structural properties 
of highly nonlinear collapsed objects (such as galaxies and galaxy 
clusters), have been widely investigated in the recent past by several 
works mostly based on dedicated large N-body simulations (see 
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e.g. Maccio et al. 2004; Baldi et al. 2010; Li & Barrow 2011 ; Baldi 
2011a, 2012b; Cervantes et al. 2012; Marulli, Baldi & Moscardini 
2012; Giocoli et al. 2013; Moresco et al. 2014; Carlesi et al. 
2014a,b). However, in these rich and high-dense environments the 
effects produced by cDE are expected to be significantly modi¬ 
fied by the complex and not yet fully understood baryonic pro¬ 
cesses occurring within astrophysical objects. Therefore, studying 
the properties of underdense regions of the universe might repre¬ 
sent a complementary approach to investigate cDE scenarios, and 
might provide a direct test of such cosmological models through a 
direct comparison with the properties of the observed cosmic voids 
(hereafter CV). 

Although the existence of CV - defined as large underdense 
regions of the Universe - was one of the earliest predictions of the 
standard cosmological model (Hausman, Olson & Roth 1983), and 
the observational discovery of CV dates back to over than 30 years 
ago (see Gregory, Thompson & Tifft 1978; Kirshner et al. 1981), 
it is only in recent years that systematic studies about CV have be¬ 
come possible thanks to the increasing depth and volume of current 
galaxy surveys and to the advent of large numerical simulations that 
allow to predict with high accuracy the topology of the cosmic web. 

The recent interest for CV is mostly related to their yet 
unexploited potential to probe cosmological models and con¬ 
strain cosmological parameters, thanks to the claimed univer¬ 
sality of their general statistical and structural properties (see 
e.g. Colberg et al. 2005; Ricciardelli, Quilis & Planelles 2013; 
Ricciardelli, Quilis & Varela 2014). In particular, CV might repre¬ 
sent a population of ideal spheres with a homogeneous distribution 
in the Universe at different redshifts, so that their size evolution can 
be used to characterise the expansion history of the Universe by 
means of the Alcock & Paczynski (AP) test (Alcock & Paczynski 
1979), as already pointed out by recent works (Sutter et al. 2012, 
2014). 

Furthermore, CV might have an impact on the observed prop¬ 
erties of the Cosmic Microwave Background (CMB), and the per¬ 
sisting CMB anomalies like the Cold Spot could be explained as 
resulting from the Integrated Sachs-Wolfe effect over large CV, 
as suggested by different works (Rees, Sciama & Stobbs 1968; 
Szapudi et al. 2014; Kovac et al. 2014; Finelli et al. 2014). The 
next generation of large galaxies surveys such as the ESA Euclid 
mission (Laureijs et al. 2011; Amendola et al. 2013) are expected 
to detect gravitational lensing from medium size CV with which 
it will be possible to directly constrain the void density profiles 
without resorting on luminous tracers like galaxies, which would 
require to model their bias (Izumi et al. 2013; Krause et al. 2013; 
Melchior et al. 2014; Clampitt & Jain 2014). 

CV are therefore one of the most appealing and promising cos¬ 
mological probes: being almost empty, their growth during the cos¬ 
mic history should be at most weakly nonlinear and their properties 
could be possibly affected by the nature of DE and by the properties 
of the primordial density field in which they evolve (Odrzywolek 
2009; D’Amico et al. 2011; Bos et al. 2012; Gibbons et al. 2014). 
In particular, the shape of CV has been shown to be very sensitive 
to the equation of state of the DE component (Lavaux & Wandelt 
2010). Defining the properties of CV in different cosmological 
models can then represent an important handle to discriminate be¬ 
tween these models. 

The present work focuses on the investigation of the prop¬ 
erties of CV in the standard ACDM cosmology, as well as in a 
series of competing cDE models. This has been done by extract¬ 
ing the population of CV from both the cold dark matter (CDM) 
and the halo distributions arising in large cosmological N-body 


simulations of these different cosmological scenarios. To this end, 
we made use of the publicly available data of the CoDECS sim¬ 
ulations (Baldi 2012c), including three different models of DE 
interaction besides a ACDM reference run. We identified CV in 
the Codecs runs with VIDE (Void IDentification and Examination 
toolkit, Sutter et al. 2015c), a substantially modified version of the 
publicly-available void finder ZOBOV (ZOnes Bordering On Void¬ 
ness, Neyrinck 2008), and compared the statistical and structural 
properties of the resulting void catalogs. Our results show that cDE 
models are characterised by an excess of large CV in the CDM 
distribution with respect to the reference ACDM cosmology, as ex¬ 
pected from their higher normalisation of linear perturbations at 
low redshifts. This is consistent with the theoretical predictions 
on the abundance of CV presented in Pisani et al. (2015), while 
the latter work seems to be in contrast with the recent findings of 
Sutter et al. (2015b) for the case of coupled dark energy simulations 
normalised to the same perturbations amplitude. 

Nonetheless, we also found that the differences in the CV 
properties among these different models significantly change when 
the CV are identified in the distribution of collapsed halos rather 
than in the CDM distribution itself. These results suggest that, con¬ 
trary to what has been claimed in some other recent works (see 
e.g. Sutter et al. 2015a), the bias of the tracers of the density field 
employed to identify the CV might have a significant impact on 
the possibility of using the obtained CV catalogs to probe cosmol¬ 
ogy. Therefore, in the present work we will show that a random 
subsampling of a simulated CDM distribution to match the density 
of tracers expected for any given galaxy survey does not actually 
provide a faithful representation of the discriminating power of the 
survey with respect to different competing cosmological models. 

The paper is organised as follows. In Section 2 we briefly de¬ 
scribe the cDE models considered in the present work and we recall 
the main features of the CoDECS runs. In Section 3 we describe the 
void finder algorithm and our method of analysis, and in Section 4 
we present the properties of the CV in the CoDECS simulations. Our 
conclusions are summarised in Section 5. 

2 COUPLED DARK ENERGY COSMOLOGIES 
2.1 The models 

In this work we aim at studying the statistical and structural prop¬ 
erties of CV in the context of coupled dark energy (cDE) cos¬ 
mologies. In these models, dark energy is represented by a classi¬ 
cal scalar field (j> moving in a self-interaction potential V{<j)) and 
directly interacting with CDM particles through an exchange of 
energy-momentum, quantified by a coupling function Here 

we will give only a very essential summary of the main features of 
cDE models, and we refer the reader to Amendola (2000); Baldi 
(2011b, 2012b) for a more thorough discussion. 

The background dynamics of cDE cosmologies is described 
by the set of equations: 


pf pf — 

0, 

(1) 

ph + '^Hph = 

0, 

(2) 

Pc “h SHPc = 


(3) 



(4) 


where the subscripts r, b, c and (/>, indicate the energy densities p of 
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radiation, baryons, CDM, and the dark energy field cj), respectively, 
and where the Huhble function is given as usual by 

--2 SttG , , 

H = g (Pt + Pc + pb + p<^) , (5) 

with Mpi = I/SttG being the reduced Planck mass. In the above 
equations the field 0 is expressed in units of Mpi and an overdot 
represents a derivative with respect to cosmic time while a prime 
denotes a derivative with respect to the field itself. The source terms 
on the right-hand side of Eqs. 3 and 4 define the interaction between 
the dark matter and the dark energy components, with a strength 
given by the coupling function 

At the level of linear density fluctuations, the interaction mod¬ 
ifies the gravitational instability processes that govern the evolution 
of perturbations as a consequence of a long-range ^/f/r force medi¬ 
ated by the dark energy field and acting between CDM fluid el¬ 
ements. In the Newtonian limit and on sub-horizon scales, these 
effects turn into the following set of modified linear equations 
(Amendola 2004; Pettorino & Baccigalupi 2008; Baldi 2012a): 


5c 


5b 


-2H 



5c -f 47rG[pb5b -P Pc5crc] , 


— 2H6h + 47rG[pb5b -p pc5c], 


( 6 ) 

(V) 


where pi represents the background density of the i-th fluid and 
5i = Spifpi its density perturbation. The factor Tc = 1 -P 4/3c/3 
represents the additional fifth force appearing only in the CDM 
equation while the term /3c0 is a velocity-dependent acceleration 
arising as a consequence of momentum conservation. Similar ad¬ 
ditional terms characterise the interaction among a discrete set of 
CDM particles in the non-linear regime (see Baldi et al. 2010). 


2.2 The Codecs simulations 

For our investigation we will make use of the publicly available 
data of the CoDECS simulations (Baldi 2012b) that represent the 
largest suite of cDE simulations to date. These simulations are car¬ 
ried out with a suitably modified version of the TreePM N-body 
code GADGET (Springel 2005) that self-consistently implements 
all the above mentioned effects characterising cDE cosmologies 
(Baldi et al. 2010). 

For the present work, we will employ the outputs of the L- 
CoDECS simulations, which follow the evolution of 1024® CDM 
particles and as many baryonic particles in a periodic cosmological 
box of 1 comoving Gpc/h a side. Both CDM and baryonic particles 
are treated as collisionless particles, but they experience different 
accelerations as a consequence of the interaction between the CDM 
and the dark energy fields. 

The Codecs suite includes six different cosmological models, 
for of them are considered in this paper: the reference ACDM cos¬ 
mology, a cDE model (EXP003) characterised by a constant posi¬ 
tive coupling Pc > 0 and an exponential self-interaction potential 
of the form V{(j)) = A exp (—ap), a further model (EXP008e3) 
with the same potential but with an exponential coupling, PdP) = 
Pq exp (Pip), and a final scenario (SUGRA003) with a constant 
negative coupling, pc < 0 and a SUGRA (Brax & Martin 1999) 
self-interaction potential Vlp) — exp (—0^/2). A sum¬ 

mary of the models parameters is shown in Table 1. All the models 
have the same amplitude of perturbations alz = zcmb , resulting in 
a different amplitude of linear density perturbations at the present 
epoch (and consequently different values of as). 

In the present work we will also make use of the public 


halo catalogs of the CoDECS simulations, that have been generated 
through a Friend-of-Friend (FoF) algorithm with a linking length 
of 0.2 times the mean inter-particle separation. 


3 METHODOLOGY 

We employ the publicly available void finder VIDE (Sutter et al. 
2015c) to identify CV in the CDM and halo distributions extracted 
from the snapshots of the CoDECS simulations within the different 
cosmological models described above. VIDE embeds the ZOBOV 
algorithm, which allows to identify depressions in the density dis¬ 
tribution of a set of points. In the following, we provide a very 
short summary of how ZOBOV works, and we refer to the original 
ZOBOV paper (Neyrinck 2008) for a more detailed discussion. 

Firstly, ZOBOV associates a cell to each tracer (a CDM par¬ 
ticle or a halo) using a Voronoi tessellation scheme, i.e. the cell 
c associated to the particle (or halo) p is defined as the region 
of the box which is closer to p than to any other particle (or 
halo) in the box. Secondly, the algorithm identifies local density 
minima among these cells: a density minimum is defined as a 
Voronoi cell with a lower density (i.e. a larger volume) than all 
other cells around it. Thirdly, ZOBOV joins together the Voronoi 
cells surrounding a local density minimum until cells with larger 
and larger density are found, and it identifies CV as the union 
of these cells. CV are joined together via the Watershed Trans¬ 
form (see Platen, van de Weygaert & Jones 2007), which naturally 
creates a hierarchy in the structures of CV. All these procedures 
are performed also by the ZOBOV version included in the VIDE 
toolkit. Additionally, VIDE provides several different void catalogs 
for which various types of sample selections (as e.g. different cuts 
on the void density contrast or on the void central overdensity) are 
applied on top of the original ZOBOV sample. In particular, as CV 
are found to define a complex hierarchy, with smaller voids being 
embedded in larger ones, VIDE provides for each identified void 
the corresponding hierarchy level, and according to this classifica¬ 
tion a sample of main CV (i.e. those CV that are not embedded 
in larger voids and that represent the top of their own void hierar¬ 
chy) is produced. We employed a slightly modified version of this 
selection procedure (see Baldi et al. in prep.) to remove pathologi¬ 
cal CV from the catalog and obtain a more statistically robust and 
convergent sample of main CV. 

Finally, since local density minima can also be found in 
overdense regions, we decide to remove from the main void 
catalogs the CV with a density minimum larger than 20% of 
the mean density of the Universe (which is one of the standard 
cuts provided by VIDE), in order to select only well defined 
CV for our comparisons. For each identified void, ZOBOV also 
calculates the probability that the void might arise in a uniform 
Poissonian distribution of points, which is directly related to the 
density contrast between the minimum density of the void and its 
boundary. As this density contrast is provided for each void also 
by the VIDE catalog, we remove CV with a density contrast below 
1.57, corresponding to a probability of arising as Poisson noise 
larger than 2a (see Neyrinck 2008). 

VIDE defines CV as spherical regions centered in the barycen- 
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Table 1. A summaiy of the cosmological models investigated in the present work and their main parameters. See Baldi (2012b) for details. 


Model 

Potential 

a 


II 

0 - 8(2 = 0) 

ACDM 

VW = A 

- 


-1.0 

0.809 

EXP003 

V{4>) = Ae-°‘^ 

0.08 

0.15 

-0.992 

0.967 

EXP008e3 

V{<l>) = 

0.08 

0.4exp[3<()] 

-0.982 

0.895 

SUGRA003 

V{(t>) = 

2.15 

-0.15 

-0.901 

0.806 


ter, Xc, of the underdense regions provided by ZOBOV, where: 

N 

X. = -, ( 8 ) 

T.vr 

i=l 

and xf and Vi are the positions of the i — th tracer and the volume 
of its associated Voronoi cell, while N is the number of tracers 
included in the void. The radius of the sphere (i.e. the effective 
radius of the void, Res) is then computed from the overall volume 
of the underdense region by assuming sphericity: 

N 

VvoiD ■ (9) 

1=1 

It has been shown that different void finders based on dynam¬ 
ical criteria, instead of density or geometry criteria, might reduce 
the shot noise error (Elyiv et al. 2015). Nevertheless, the CV finder 
used here is accurate enough for the purpose of the present anal¬ 
ysis, as we investigate the main properties of large CV extracted 
from dense numerical simulations. 


4 THE STATISTICS OF VOIDS IN THE CODECS 

With the catalogs of CV extracted from the CoDECS simulations 
as described in the previous Section at hand, we perform some ba¬ 
sic analyses of the statistical and structural properties of the CV 
in the different cosmologies, namely the CV filling factor (i.e. the 
fraction of the cosmic volume occupied by CV), their size distri¬ 
bution (i.e. the abundance of CV as a function of their size), and 
their stacked radial density profiles, and compare these observables 
to the reference ACDM case (see also Li 2011). We perform such 
comparison for CV identified both in a randomly subsampled CDM 
density field and in the distribution of collapsed halos, to highlight 
how the use of tracers with different bias might result in a different 
relative behaviour of the models. 

4.1 Void statistics in the CDM distribution 

Let us start by considering the CV catalogs extracted from the 
CDM density field, i.e. directly from the CoDECS snapshots at dif¬ 
ferent redshifts. To better handle the simulation data we have made 
use of the subsampling routine included in VIDE to randomly sub¬ 
sample the CDM particles of the simulation snapshots down to 
an average density of 2 x particles per cubic h~^ Gpc. For 
this comparison we will focus only on two out of the four models, 
namely the reference ACDM cosmology and the EXP003 scenario, 
which is the most extreme realisation (in terms of deviations at the 
background and linear perturbations level) of cDE models that we 
have at our disposal. 

First of all, we compare the evolution of the volume fraction 
occupied by CV at different cosmic times, also known as the void 
filling factor, to check whether the interaction between DE and 



Figure 1. The redshift evolution of the volume fraction of CV identified in 
the CDM distribution for the ACDM (black solid line) and EXP003 (blue 
dashed line) models. The shaded area (shown only for ACDM) represents 
the uncertainty, computed with the jackknife method. 


CDM particles implemented in our extreme cDE model has an im¬ 
pact on such fraction. Fig. 1 displays the evolution of the CV vol¬ 
ume fraction, Vvoids/Vtot. where VVoids is the sum of all the 
main CV volumes in a given snapshot of the simulation and Vtot 
is the total volume of the box, i.e. 1 h~^ Gpc®. The statistical error, 
shown in Fig. 1 by the shaded grey region around the ACDM line, 
has been computed with a jackknife method. 

As expected, the volume fraction of CV increases with time, 
irrespectively of the underlying cosmological model, due to gravi¬ 
tational instability. Moreover, as one can see from Fig. 1, the vol¬ 
ume fraction occupied by CV in the cDE model EXP003 is signif¬ 
icantly larger than in the reference ACDM cosmology, reflecting 
the higher normalisation of the amplitude of linear perturbations at 
low redshifts in EXP003. More quantitatively, the volume fraction 
in EXP003 is roughly 40% larger than the corresponding ACDM 
fraction, at all redshifts between 2 = 1 and 2 = 0. Clearly, the ob¬ 
served differences between the cDE model and the standard ACDM 
cosmology are statistically significant. 

As a second step, we compare the relative abundance of CV 
as a function of their size, by computing in the two cosmological 
models the differential size distribution (hereafter DSD), defined 
as the number of CV with an effective radius J?eff falling within a 
set of size bins. In the upper panels of Fig. 2 we show the DSD at 
three different redshifts (2 = {0,0.55,1}, from left to right) for 
the two cosmological models (black solid lines for ACDM and blue 
dashed lines for EXP003), while in the bottom panels we show the 
relative difference with respect to the reference ACDM cosmology 
(in units of its statistical error a). As one can see in the figure, at 
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2 = 0 the number of CV in the cDE cosmology with i?eff ^ 20 
h~^ Mpc is at least 50% larger than in ACDM; this difference cor¬ 
responds to more than 4 (t. At 2 = 0.55 the same ratio applies to 
CV with i?eff > 15 h~^ Mpc, and at 2 = 1 to CV with R^s > 12 
h~^ Mpc, with differences corresponding to 5cr and 7a, respec¬ 
tively. Therefore, also in this different statistics the two models are 
clearly distinguishable from each other. 

As a third statistic of our CV samples, we investigate the 
average stacked density profiles of CV having a comparable size. 
Many recent works (e.g. Ricciardelli, Quilis & Planelles 2013; 
Ricciardelli, Quilis & Varela 2014; Hamaus, Sutter & Wandelt 
2014) suggested that the average profile of CV is self-similar in 
the standard ACDM cosmology, which makes CV an ideal target 
for geometrical tests such as the AP test. Therefore, we now aim 
at investigating whether the interaction between DE and CDM 
particles might induce some additional features on the density 
profile of CV. 

To this end, we first compute the spherically-averaged radial 
density profile of each individual void by estimating the CDM 
density within a series of logarithmically equispaced spherical 
shells centred in the barycentre of each void and normalised to the 
void effective radius 7?eff. The profiles are then stacked for CV 
with similar i?eff. Since the profiles of each void is calculated in 
units of Tieff in the first place, the stacking procedure basically 
consists in the calculation of the mean profile in each logarithmic 
radial bin. We randomly included 100 CV for each bin. The 
results are presented in Fig. 3, where we show the comparison 
of the stacked density profiles obtained using the CDM void 
catalog in the standard ACDM cosmology and in the EXP003 
model. In the upper panels the error bars represent the corrected 
sample standard deviation computed on the 100 randomly se¬ 
lected CV. The shape of the profiles is qualitatively the same as 
found in previous works (Ricciardelli, Quilis & Planelles 2013; 
Ricciardelli, Quilis & Varela 2014; Hamaus, Sutter & Wandelt 
2014): we observe a deep underdensity at i? —>■ 0 and a compen¬ 
sative overdensity at i? —>■ Res- At i? > 1.5 • ReS the profiles 
reach the mean density of the Universe. In the lower panels, 
we plot the relative difference between the models in units of 
the statistical significance a computed as the sample standard 
deviation propagated to the relative difference. The grey shaded 
area represents a Tier significance. 

As the figure clearly shows, the stacked profiles of EXP003 do 
not show significative differences from ACDM at the considered 
redshifts, with deviations always lying well within the mean square 
en'or. Nonetheless, we can observe that, at ii —>■ 0, EXP003 gener¬ 
ally shows a density 10 — 25% smaller than ACDM, thereby show¬ 
ing that CV are emptier in cDE. On the other hand, the compen¬ 
sative over-density at i? = i?eft for EXP003 looks generally more 
prominent than ACDM (except for CV with 5 < Rea\h~^ Mpc ] < 
8). Therefore, although with a low statistical significance, the cen¬ 
tral regions of the main CV appear to be more underdense in cDE 
models than in ACDM, which is expected to result in a correspond¬ 
ing stronger signal in void lensing surveys. 

4.2 Void statistics in the halo distrihution 

We will now repeat the same three statistics of the CV properties 
discussed in Section 4.1 for the CV catalogs obtained by running 
VIDE on the distribution of FoF halos extracted from the CoDECS 
simulations at the same three redshifts investigated before (i.e. 
2 = {0,0.55,1}). The use of the FoF halos as tracers of the matter 
distribution has the appealing property to mimic real observations, 


where CV are identified in the distribution of luminous galaxies. 
In particular, we have made use of the publicly available CoDECS 
halo catalogs that have been obtained through a FoF algorithm with 
a linking length 0.2 times the mean inter-particle separation. As we 
will show below, the differences between the cDE model EXP003 
and the standard ACDM cosmology in all the three statistics are 
much weaker than what previously found for the CDM distribu¬ 
tion. For this reason, we will include in this comparison also other 
two cDE models available within the CoDECS suite, namely the 
EXP008e3 and the SUGRA003 models (introduced in Section 2.2), 
in order to verify whether different realisations of the cDE scenario 
might have a stronger impact on the CV defined by the FoF halo 
distribution than the EXP003 model. Our comparison will show 
that this is actually not the case, as expected from the fact that 
EXP003 is the most extreme of the CoDECS models in terms of 
background and linear deviations from ACDM. 

First of all, in Fig. 4 we compare the CV filling factor for 
these new void catalogs, as already done in Fig. 1 for the CV in the 
CDM distribution. The void volume and the dispersion indicated 
by the grey shaded area are computed as outlined above. The figure 
shows, as expected, that the void volume fraction increases in time, 
and that the ACDM model has generally the lowest volume fraction 
with respect to the other cDE models that are characterised by a 
higher normalisation of the linear power spectrum. Nonetheless, as 
the figure clearly shows, these differences are now much smaller 
and lie within the 3 ct statistical dispersion so that no significant 
differences in the CV filling factor appear among the various cDE 
models and the standard ACDM cosmology at all redshifts. This 
result is starkly different from what found for the CV in the CDM 
distribution for the EXP003 model. 

In Fig. 5 we then display the DSD for CV identified in the 
distribution of FoF halos, analogously to what done in Fig. 2 for 
the CDM distribution. The DSD is shown for the different mod¬ 
els (ACDM by black solid line, EXP003 by blue dashed line, 
EXP008e3 by orange dot-dashed line and SUGRA003 by red dot¬ 
ted line) in the upper panels, while the bottom panels report the 
percent deviation in units of the statistical significance cr from the 
reference ACDM case. While at 2 = 0 and 2 = 0.55 no significant 
differences appear among the models, one can observe an excess 
of small CV for the EXP003 model at 2 = 1. At this redshifts 
EXP003 shows ~ 50% more CV with 7?eff < 30 Mpc than 
ACDM, although within a confidence of Icr. 

This is again a very different result with respect to what previ¬ 
ously found for the CV identified in the CDM distribution, where 
the largest differences with respect to the standard cosmological 
model appeared at the large size tail of the distribution. It should 
however be noticed that due to the different density of the tracers 
between the subsampled CDM distribution adopted in the previ¬ 
ous Section and the FoF halo distribution shown here, the mean 
separation between particles and hence the average size of CV is 
different in the two cases. Therefore, the size range that appeared 
as the large-size tail for the CDM CV (20 < i?eff [h~^ Mpc ] < 30) 
is now representing the small-size part of the void samples of the 
FoF halo distribution. For this reason, the two results might still 
appear consistent with each other despite their different qualitative 
trends. To address this issue, in Section 4.3 below we will compare 
the DSD of CV identified in a different random subsample of the 
CDM particles distribution with the same density of tracers as the 
FoF halo catalogue. Nonetheless, the clear differences between the 
background evolution of ACDM and cDE models (see Fig. 2) are 
not expected to be detected and do not appear for CV in halos. 

Before moving to this additional comparison, we conclude 
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Figure 2. Top panels: the size distribution of CV in the CDM distribution for the ACDM (black solid line) and EXP003 (blue dashed lines) models. Bottom 
panels: the relative differences between the two models in units of the standard deviation a, computed for the ACDM model. 




Figure 3. The stacked profiles of CV in the CDM distribution for the ACDM (black solid lines) and EXP003 (blue dashed lines) models. Results are displayed 
at two different redshifts, z = 0 and 2 = 1 (left and right blocks of panels, respectively) for four ranges of labeled. The error bars indicate 

the coiTected sample standard deviation in each radial bin computed on the 100 randomly selected CV, while the sub-panels display the relative difference 
between the profiles in units of the statistical significance of the averaged profile. 


our analysis of the CV extracted from the FoF halo distribution 
by comparing the void stacked density profiles as we did in Fig. 3 
for the CV in the CDM distribution. In Fig. 6 we show the anal¬ 
ogous to Fig. 3 for these new CV samples at 2 = 0 and 2 = 1. 
For both redshifts we do observe significative deviations from the 
ACDM profile only in the inner part of the CV, where shot noise 
as well as an offset of the computed centre of the void from the 
real centre can strongly affect the profiles, as pointed out by e.g. 
Nadathur et al. (2015); Nadathur & Hotchkiss (2015). We also ob¬ 
serve that the over-compensative region around ~ is not as 


prominent as in the CDM CV (see Fig. 3), once again showing 
differences between tracers of density. We notice that the density 
minimum does not lie exactly at i? ~ 0, which was not the case 
for CV in CDM. These last problem might be caused by the defini¬ 
tion of centre as in eq. 8, as pointed out by Nadathur et al. (2015); 
Nadathur & Hotchkiss (2015). Eq. 8, anyway, is sufficient for our 
purpose because, by using it, we can easily compare our results 
with other works, and because we are not discussing features re¬ 
lated to the correlation function of CV (in which the position of the 
centre of voids plays a fundamental role). 
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Figure 4. The redshift evolution of the volume fraction of CV identihed 
in the halo distribution for different cosmological models: ACDM (black 
solid line), EXP003 (blue dashed line), EXP008e3 (dot-dashed orange line), 
SUGRA003 (red dotted line). The shaded area (shown only for ACDM) 
represents the uncertainty, computed with the jackknife method. 


4.3 The impact of halo bias 

As introduced in Section 4.2, it is interesting to compare the abun¬ 
dance of CV as a function of their effective radius, 7?eff, for CV 
samples extracted from the FoF halo catalogs and from a random 
subsampling of the CDM distribution having the same number of 
tracers in the simulation box as the number of FoF halos. This will 
ensure that the mean inter-particle separation of the two samples 
of tracers is the same such that the corresponding average size of 
CV will be comparable in the two cases, thereby allowing for a 
direct comparison of the statistical properties of the two CV cata¬ 
logs over a similar range of void sizes. This approach has been fol¬ 
lowed in several recent works (see e.g. Hamaus, Sutter & Wandelt 
2014; Sutter etal. 2015b; Pisani et al. 2015), where the discrimi¬ 
nating power of CV in future galaxy surveys has been inferred from 
the expected properties of CV identified in a random subsample of a 
simulated CDM distribution having the same density of the survey 
under investigation. To this end, we have randomly subsampled the 
CDM distribution of the CoDECS snapshots at the relevant redshifts 
to a total number of 1.5 • 10® particles, corresponding to the total 
number of objects in the FoF halo catalog of the ACDM simulation 
at « = 0. 

In Fig. 7 we show the equivalent to Fig. 2 for this new random 
subsampling and compare the abundance of CV in the ACDM and 
EXP003 models. We observe that, although the range of void sizes 
is now comparable to what was shown in Fig. 5, the comparison 
between ACDM and EXP003 still appears starkly different in the 
two cases. Also in this case, as already shown for a denser sam¬ 
ple of CDM tracers, the EXP003 scenario includes a larger number 
of CV of large sizes at all redshifts with respect to the standard 
ACDM cosmology, with a qualitatively different trend with respect 
to what shown in Eig. 5. The comparison of Eigs. 5 and 7 clearly 
indicates that CV in the CDM distribution and CV in the distribu¬ 
tion of halos are characterised by different statistical properties. As 
the density of the two tracers is the same, these different properties 
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0.55 

1.584 
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0.733 
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1.00 

2.049 

0.504 

1.033 

1.633 

0.595 

0.972 

1.60 

2.903 

0.398 

1.155 

2.235 

0.468 

1.046 

2.00 

3.630 

0.348 

1.263 

2.739 

0.408 

1.118 


Table 2. The bias b(z) and the normalisation of the linear perturbations 
amplitude o's(z) for the ACDM and EXP003 cosmologies. The rightmost 
column for each model displays the combination b{z) ■ as{z), showing how 
this combination is much similar for the two models as compared to 0-3 
alone. As a consequence, the differences in the void populations extracted 
from the biased tracers within the two scenarios are significantly suppressed 
with respect to the case of the CV in the CDM distribution. 

must be associated with the different bias of the two samples with 
respect to the underlying true density field: while a random subsam¬ 
pling of the CDM distribution is an unbiased tracer of the density 
field, halos are biased and the bias is expected to evolve differ¬ 
ently in cDE models than in ACDM (Marulli, Baldi & Moscardini 
2012; Moresco et al. 2014). More quantitatively, the lower bias of 
the EXP003 model compensates for the higher value of the pertur¬ 
bations amplitude. In Table 2 we display the value of the bias (as 
computed in Marulli, Baldi & Moscardini 2012) and of as at vari¬ 
ous redshifts for the two models. As one can see from the last col¬ 
umn, the combination b{z) ■ as{z) is substantially closer between 
the two models compared to the value of as alone. This result sug¬ 
gests that the assumption (implicitly adopted in many recent works) 
that the properties of CV in a subsampled set of CDM particles ex¬ 
tracted from a cosmological simulation can faithfully reproduce the 
statistics of CV identified in a galaxy survey is not valid. 


In order to further validate this result, we compute the DSD 
of CV identified in the distribution of EoE halos with masses 
M > 5-10^^ Mq , thus considering tracers with larger masses and, 
therefore, with higher bias. This comparison is shown in the upper 
panels of Eig. 8, while the bottom panels display the deviation (in 
units of a) between the models. At 2 = 0 we do not observe any 
significative difference between models (in agreement with Eig. 5), 
while at larger redshifts we find that the EXP003 cDE model fea¬ 
tures a larger number of small CV (30 < Reff[h~^ Mpc ]< 60) as 
compared to ACDM, though the effect is small. The comparison 
between Pigs. 2 and 8 indicates again that CV in the CDM distri¬ 
bution and in the distribution of halos are characterised by different 
statistical properties. This result clearly shows that the bias of the 
tracers from which CV are identified has a non-trivial impact on the 
relative statistical properties of the CV sample between two com¬ 
peting cosmological scenarios. Therefore, when comparing CV in 
ACDM and cDE models, CV in halos (which are biased tracers of 
the underlying density field) will not provide a faithful representa¬ 
tion of how the models might differ in the properties of CV in the 
CDM distribution. 


5 DISCUSSION AND CONCLUSIONS 

In this paper we analysed the statistical properties of CV in ACDM 
and cDE models. In particular, we compared the properties of CV 
detected in the distribution of CDM and in collapsed halos, by 
means of a suite of large cosmological simulations, the CoDECS. 
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Figure 5. Top panels: the size distribution of CV in the halo distribution for the ACDM (black solid line), EXP003 (blue dashed lines), EXP008e3 (orange 
dot-dashed lines) and SUGRA003 (red dotted lines) models. Bottom panels: the relative differences between the cDE models and the ACDM one, in units of 
the standard deviation cr, computed for the ACDM model. 




Figure 6. The stacked profiles of CV in the halo distribution for the ACDM (black solid lines), EXP003 (blue dashed lines), EXP008e3 (orange dot-dashed 
lines) and SUGRA003 (red dotted lines) models. Results are displayed at two different redshifts, z = 0 and 2 = 1 (left and right blocks of panels, respectively) 
for four ranges of iJeff > ^s labeled. The eiTor bars in the upper panels are computed as for Fig. 3, and the sub-panels display again the relative difference of the 
profiles with respect to the ACDM one in units of the statistical significance of the averaged profile. . 


We focused on three CV statistics: the filling factor, the size distri¬ 
bution and the stacked density profiles. 

In Section 4.1 we investigate the properties of CV in the CDM 
distribution, considering the ACDM and the cDE model EXP003, 
which represents - among the available CoDECS models - the most 
extreme case showing the largest discrepancies with respect to 
ACDM in several other observables (see e.g. the results of Baldi 
2012b; Lee & Baldi 2011; Marulli, Baldi & Moscardini 2012; 
Beynon et al. 2012; Cui, Baldi & Borgani 2012; Giocoli et al. 


2013; Carbone et al. 2013; Moresco et al. 2014; Pace et al. 2015; 
Giocoli etal. 2015). Our main results can be summarised as fol¬ 
lows. 

(i) The filling factor of CV detected in the CDM distribution in 
the EXP003 model is significantly larger than in the ACDM case, 
as expected due to the higher normalization of the amplitude of 
linear perturbations at low redshift (Fig. 1). More quantitatively, the 
volume fraction in EXP003 is ~ 40% larger than the corresponding 
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Figure 7. Top panels: the size distribution of CV in the CDM distribution diluted to the same density of the halo catalogue, for the ACDM (black solid lines) 
and EXP003 (blue dashed lines) models. Bottom panels: the relative differences between the two models, in units of the standard deviation, cr, computed for 
the ACDM model. 



Figure 8 . Top panels: the size distribution of CV identified in the distribution of halos with mass M > 5 ■ 10^^ Mq, for the ACDM (black solid lines), 
EXP003 (blue dashed lines), EXP008e3 (orange dot-dashed lines) and SUGRA003 (red dotted lines) models. Bottom panels: the relative differences between 
the cDE models and the ACDM one, in units of the standard deviation cr, computed for the ACDM model. 


ACDM fraction: based on a jackknife approach, this is detectable 
with a very high statistical significance. 

(ii) For what concerns the DSD (Fig. 2), we found an ex¬ 
cess of large CV in the EXP003 model with respect to the refer¬ 
ence ACDM cosmology, consistently with the general findings of 
Pisani et al. (2015). Quantitatively, the excess is around 50% with 


a difference larger than 4a. The radius at which this excess starts 
to be significant decreases with redshift, being R ~ 20,15,12 
h~^ Mpc at z = 0, 0.55, 1, respectively. 

(iii) The shape of the stacked density profile (Fig. 3) is qualita¬ 
tively similar to what previously found in the literature (i.e density 
minima around the centres of the CV and overdense compensation 
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regions at i? ~ Re«)- The void profiles in cDE models are not 
significantly different from what observed in the standard cosmol¬ 
ogy. Nonetheless, we can observe that close to the CV centres, the 
EXP003 model generally displays a density 10 — 25% smaller than 
ACDM, thus showing that CV tend to be emptier in cDE models. 
On the other hand, the compensative overdensity at i? ~ ReS in the 
EXP003 case looks more prominent than in ACDM. All of these 
features are expected considering that the evolution of the back¬ 
ground perturbations in cDE scenarios is faster than ACDM due to 
the fifth force associated with the coupling. 

In Section 4.2 we then focused on CV identified in the halo 
distribution, finding that the comparison between cDE models and 
the reference cosmology is very different from what found for 
CDM. More specifically, we find the following results: 

(i) The filling factor of CV in halos is not strongly dependent 
on the considered cDE model (Pig. 4). Only minor, not significant 
differences are found in the volume fractions at all redshifts con¬ 
sidered in this paper. This last result is starkly different from what 
observed in CV samples detected in the CDM distribution. 

(ii) The comparison of DSD (Pig. 5) in the halo distribution does 

not reveal sensible differences between cDE models and the ref¬ 
erence one. This result is again substantially different from what 
found in Section 4.1 . We connect this discrepancy with the impact 
of the halo bias on CV properties. To test such effect, we compare 
DSD of CV in a random subsample of the CDM distribution with 
the same density of tracers as the FoF halo catalog. Again, in this 
last case CV in CDM do not show the same relative trend in the 
DSD as for the CV in the halo distribution (Fig. 7). The impact of 
the halo bias can be observed also by increasing the minimum mass 
of halos used as tracers: in Fig. 8 we show that including only halos 
with large masses (> 5 • 10^^ • Mq) the cDE models show an ex¬ 
cess of CV with 30 < Mpc ]< 60 at z = 0.55,1, which 

is not seen in Fig. 5. 

(iii) The density profile of CV in halos does not look like an ef¬ 
fective probe to discriminate among cDE models. Indeed, as shown 
in Fig. 6, the stacked profiles of CV in cDE models are only 
marginally distinguishable from the ACDM case, and only in the 
very innermost parts. 

To conclude, the main result of this work is that the properties 
of CV in different cosmological models are strongly affected by the 
choice of the tracers of the underlying density field used to detect 
them (halos or CDM particles): this is caused by the impact of the 
halo bias on the structural properties of CV: as the bias evolves dif¬ 
ferently for different cDE models, this is reflected in a non-trivial 
way on the properties of the associated CV sample. Our results indi¬ 
rectly challenge the assumption made in several recent works that a 
subsampled distribution of simulated CDM particles with the same 
density of the expected tracers of a real galaxy survey might pro¬ 
vide reliable predictions about the effective discriminating power 
of CV in that survey. 
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